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ABSTRACT 

Neutralino annihilation and elastic scattering cross sections are derived which 
differ m important ways from previous work. These are applied to relic abun- 
dance calculations and to direct detection of neutralino dark matter from the 
galactic halo. Assuming the neutralino is the lightest supersymmetric particle 
and that it is less massive than the Z°, we find relic densities of neutralinos 
greater than 4% of critical density for almost all values of the supersymmetric 
parameters. We constrain the parameter space by using results from PETRA 
(chargino mass less than 23 GeV) and ASP, and then assuming a critical density 
of neutralinos, display event rates in a cryogenic detector for a variety of models. 
A new term implies “spin independent” elastic scattering even for these majo- 
rana particles and inclusion of propagator momenta increases detection rates by 
10% to 300% even for pure photinos. Z°-squark interference leads to very low 
detection rates for some values of the parameters. The new term in the elastic 
cross section dominates for heavy, mostly spinless materials and mitigates the 
negative interference cancellations in light materials; except for the pure photino 
or pure higgsinos cases where it does not contribute. In general, the rates can 

be substantially different from the pure photino and pure higgsino special cases 
usually considered. 
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I. Introduction 


One of the most interesting unsolved problems in physics today is the iden- 
tity of the dark matter (DM) known to exist in galactic halos. Several lines of 
reasoning lead one to hypothesis that the DM may not be ordinary baryonic 
material but rather may consist of some, as yet undiscovered, elementary par- 
ticle. If so, it is likely that substantial quantities of these particles exist in our 
galaxy’s halo at the position of the Earth and it may be possible to detect these 
particles as they pass through detectors in laboratories on Earth. Indirect 
methods, such as detecting the products of DM particle-antiparticle annihilation 
have also been suggested, the most promising of which use the enhancement in 
density which results from capture of DM particles into the body of the Sun or 
Earth . 5,6 

Crucial to all detection schemes are the DM particle cross sections. A zero 
elastic scattering cross section off ordinary matter would mean no possibility 
of direct direction and no capture of particles into the Earth or Sun. Particle- 
antiparticle annihilation cross sections determine, in part, the relic abundance 
of DM particles and also the flux of detectable annihilation products. Since 
there is little astrophysical data which bear on these cross sections we must turn 
to particle physics models to make predictions. A wide spectrum of candidate 
particles have been proposed, but the most interesting are those which were 
originally proposed for non-astrophysical reasons, and only subsequently turned 
out to be suitable as DM candidates. The most popular particles in this class 
are the neutrino, the axion, and the lightest supersymmetric particle (LSP). In 
this Report we consider only the LSP. 

Supersymmetry has been immensely popular among theorists for the past 
decade. It seems to be an essential ingredient of theories which unify gravity 
with the other interactions and in its low-energy manifestation provides an ele- 
gant solution to the hierarchy problem. Supersymmetry has also been popular 
among experimentalists, but although much effort has been expended in the 
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search, no supersymmetric particles have been discovered. A major goal of the 
next generation of particle accelerators will be to discover or set limits on super- 
symmetric particles. In the minimal supersymmetric extension of the Standard 
Model which we consider here, there is a new fermionic (bosonic) superpartner 
for every bosonic (fermionic) ordinary particle, so there is a doubling of the num- 
ber of Standard Model particles. The masses of these particles are undetermined 
in the most general scheme, but if supersymmetry is to solve the hierarchy prob- 
lem, some superpartners must exist with masses below roughly a TeV. In most 
versions of these models there is also a conserved multiplicative quantum number 
= — 1 2S+3B+L R.panty implies that the lightest superpartner is stable 

and therefore a DM candidate. The identity and cosmology of the LSP was con- 
sidered in 1984 by Ellis et al., 9 and since then, the lightest neutralino (*), a linear 
combination of the photino, zino and two neutral higgsinos, has been considered 
the likely candidate for LSP. In considering detectability, most authors have, 
however, considered only the pure photino and pure higgsino, two special cases 
of the general neutralino. If the neutralino is very light, then one might expect 
a reasonably pure photino or higgsino, but there are no strong theoretical or ex- 
perimental reasons to expect such a light LSP, and as the mass of the neutralino 
increases, a pure photino or higgsino becomes more and more unlikely. 

One of the most remarkable characteristics of neutralinos is that over most 
of the supersymmetric parameter space a relic density p^ near critical density 
is obtained. We emphasize that a value of D* (defined as p-il Pent) near 1, is 
therefore very natural, and this fact adds motivation to the search for neutralino 
dark matter (See also Ref. 9). In fact, even if neutralinos do not constitute the 
DM responsible for galactic rotation curves, we still expect > .04 in our galac- 
tic halo as long as they exist and are the LSP. So if low-energy supersymmetry 
does exist we expect to find a substantial fraction of the mass of the universe in 
neutralinos and it is therefore worth making the effort to detect them. 

In this Report we recalculate the cross sections for neutralinos and find several 
important differences with respect to earlier work. Using these cross sections 
we show that direct detection event rates in cryogenic detectors can be very 
different from the pure photino and pure higgsinos special cases. We also find 
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that as the neutralino becomes more massive pure photinos and especially pure 
higgsinos become very unlikely. Including a new term in the elastic scattering 
cross section and using a method due to Shifman tt al. y which includes the 
effects of heavy quark loops, we find a piece of the elastic cross section which 
is proportional to the mass of the scattering nucleus. This means neutralinos 
might be detectable with mostly spinless material such as mercury, germanium 
or helium, even though majorana particles are usually thought to have only spin 
dependent” interactions. 

Inclusion of several new terms in the annihilation cross section, as well as 
a different sign interference term, changes relic abundance results for neutrali- 
nos, and the inclusion of previously ignored effects such as propagator momenta 

changes the results even for pure photinos. Finally, since the same couplings 

, ..12 

are involved in the elastic, annihilation, and production cross sections we can 

extract information from various accelerator experiments. In particular, large 

33 

areas of neutralino parameter space are ruled out by experiments such as UA1, 
PETRA, 14 and ASP, 15 and we find that neutralinos with mass less than about 
5 GeV are unlikely. 

The plan of this Report is as follows. In Sec. II we define the models un- 
der consideration and the parameter space to be explored. We then calculate 
the XX99 matrix element { q is any quark or lepton) and the annihilation cross 
section (xx gg). For use in relic abundance calculations we also present the 
non-relativistic expansion including the possibly important effect of propagator 
momenta. Exploiting crossing symmetry, we use the annihilation matrix element 
to find the elastic scattering (xq xq) and production (e+e~ -► xx) cross 
sections. For use in comparison with ASP results, we finally present the cross 
section cr(e~^e — — * XXT) the soft photon limit. 

In Sec. Ill we re-derive and extend some of the results from Sec. II using an 
effective Lagrangian approach. In particular, we find the elastic scattering cross 
section of neutralinos off nuclei rather than neutralinos off quarks and show that 
due to heavy quark loops, a “spin independent” interaction exists for the general 
neutralino. This interaction is proportional to the zino component and so does 
not contribute to pure photino or pure higgsino scattering. We also point out a 
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difference (with respect to earlier work) of the sign of an interference term. 

In Sec. IV we use the annihilation cross section to calculate the relic abun- 
dance of neutralinos and show that ft* ~ 1 is very natural. We also show that 
apart from near the Z° pole a very low relic density of neutralinos is unlikely. 
We also describe and implement several consistency requirements such as the 
neutralino being lighter than the squarks and the charginos. 

In Sec. V we use the elastic cross section to find the rate at which DM neu- 
tralinos interact with nuclei in a cryogenic detector. We show that the new 
(scalar) term can give enhancements of several orders of magnitude for very 
heavy nuclei. We also point out that with the standard “spin dependent” (axial 
vector) cross section very low event rates generically occur due to negative inter- 
ference between Z and squark exchange. These cancellations are substantially 
mitigated by the new scalar term even for light elements. ^Ve reduce parameter 
space throughout by requiring that ft* = 1. We also mention that the inclusion 
of propagator momenta can be important, especially when the neutralino mass 
nears the squark mass. This is illustrated for the pure photino case where, for a 
neutralino mass of 30 GeV , the event rate with propagator momenta included is 
a factor of two higher than without. We emphasize that the event rates for the 
general neutralino can be quite different from the pure photino and pure higgsino 
special cases, and that since detection should be aimed at the general neutralino, 
all terms in the cross section should be taken into account when selecting detector 
material. 

In Sec. VI we show how accelerator experiments can compliment dark mat- 
ter searches. In particular, the PETRA result 14 that m i± < 23 GeV eliminates 
large regions of parameter space and as these limits improve (or the chargino 
is detected) the chargino constraint will be come more and more important 
in deciding the detectability of neutralino DM. In addition, the ASP limit 15 
<r(e + e~ — > 7 -f- missing) < .03 pb is re-examined for the case where “missing” is 
a pair of neutralinos, and again parts of parameter space are eliminated. While 
fhe PETRA constraint eliminates mostly states with a substantial zino compo- 
nent, the ASP constraint rules out light, mostly photino states, and together the 
constraints rule out almost all states with mass less than about 5 GeV. 
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Finally, in Sec. VII we consider the effect of relaxing several simplifying as- 
sumptions made during the first six sections. We consider the effect the new 
EMC spin dependent structure functions 16 have on direct detection, and find 
a sensitivity to the composition of the neutralino and shell model parameteri- 
zation of the nuclei whenever the axial vector term dominates the elastic cross 
section. For shell model “neutrons” there is a general lowering of rates, except 
for photino-like states which are greatly enhanced. The change in rates for shell 
model “protons” is not as large except again for photino-like states which are 
somewhat suppressed. However, while details can vary substantially, the general 
picture remains unchanged. We also consider the effect of non-degenerate selec- 
tron and squark masses (M f and Mg respectively). We show results for Mg = 3 Af* 
which, along with the previously considered M, $ = M j-, bracket most of the mod- 
els we surveyed from the literature. Event rates are in general one or two orders 
of magnitude lower, since annihilation through the slepton channels (oc M- ) 
can be strong, resulting in weaker couplings to achieve = 1, and these weaker 
couplings remain for the elastic scattering which proceeds via the heavier squark 
(oc Mr 4 ) exchange. Sec. VIII sums up the Report. 

II. Cross sections 

In this section we calculate the matrix element and cross section for the 
annihilation (xx -* Ql) °f two neutralinos into a fermion and antifermion, and 
then use crossing symmetry to find the elastic scattering (x9 X?) an d tbe 

anomalous single photon (ASP) (e + e“ -+ XX7) cross sections. Here q stands for 
any quark or lepton and e~ is the electron. 

Throughout we will use the minimal supersymmetric extension of the Stan- 
dard Model as described in Ref. 7 (see especially the appendices) and Ref. 8. 
This is a group of models which is minimal, in the sense that it contains the min- 
imum number of new particles, and general, in the sense that it has all possible 
CP-conserving soft supersymmetry breaking terms included in the Lagrangian. 
In these models there exist four neutralinos which are linear combinations of the 
supersymmetric partners of the neutral W , the i?, and the two neutral Higgs 
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bosons. These can also be characterized as the photino, zino and two neutral 
higgsinos. Only the lightest will be stable (we assume a conserved R parity and 
also that the lightest neutralino is the LSP) and we denote it as 

X = Zn B + Z 12 W 3 + Z n H 1 + Z 14 H 2 , ( 1 ) 

where the Z\{ are the elements of the real orthogonal matrix which diagonalizes 
the neutralino mass matrix; that is, if Z\\ = Z\ 2 = 0, x is a pure higgsino, if 
Z 11 = cos 9 W , Z\2 = sin 0 W , x is a pure photino, and if Z n = - sin 0 Wi Z 12 = 
cos 0 W , x is a pure zino. The assumption of CP-conservation ensures that the Z{j 
are real. 

The neutralino masses and the Z,y’s are fully determined by four parameters: 
tan/?, fi) M , and M , where tan/? = v 2 jv\ is the ratio of Higgs vacuum expecta- 
tion values, M and M f are soft supersymmetry breaking parameters, and /z is a 
supersymmetric Higgs mass. Throughout, we make the standard 7 simplification 
M = | M tan 2 0 W which is true if the theory is eventually embedded in any sim- 
ple grand unified group. Overall then we have three undetermined parameters: 
tan/?, M } and /z, and it is this parameter space we explore. ^Ve consider several 
representative values of tan (3 and M and fi in the wide range 0 < < 1 TeV. 

For a neutralino of mass mx less than the Z° mass m^, the annihilation, 
elastic scattering, and production (e + e — > xx) processes are all given by the 
five Feynman diagrams shown in Figure 1. There are four squark (or slepton) 
exchange diagrams and one Z® exchange. The Feynman rule for the left chiral 
X 99 X vertex is — igy/ 2 (aePji + bPi) and the rule for the right chiral XQQR vertex is 
— igV 2 (aPi — €cPr) where g is the electroweak coupling constant, Pr = 1(1 +75), 
etc. and 

<Z — TTlqdq j ( ^ 

b = T 3L Zj 2 - tan 9 W (T 3L - e q )Zj 1 (2) 

c — tan 0w^qZj\ . 

Here d q = Zy 3 / cos/? for down type quarks or leptons and d q = Zj 4 / sin/? for up 
types, rriq is the mass of the quark or lepton, T3 r is the weak isospin, is the 
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charge, sin 2 9 W = . 23 , and e is the sign of the neutralino mass eigenvalue. The 
“a” factor corresponds to a higgsino type coupling, while the “6” and “c” factors 
are the photino-zino couplings. The Z°XX Feynman rule is (t'g/2 cos 0 W )(Z ^ 3 — 
Zj 4 ) 7^75 and the Z°qq rule is as usual (-ig/ cos Ow^^^lPl + crPr)- Here 
ci = T^l — e q sin 2 6 W and cr = -e q sin 2 8 W . More details on the Feynman rules 
and techniques for handling majorana fermions can be found in Ref. 7 appendices 
and Ref. 8. 

Since our results differ in the sign of interference terms and include new 
terms we present our matrix element before summing and squaring in appendix 
A. The complete summed and squared matrix element including non-degenerate 
left and right chiral squark masses is also included in appendix A. The effect of 
left-right chiral splitting has been considered previously ’ and is expected to 
be small. In order to simplify the formulas we will not consider this effect and 
set M q i = M$r = M q . In this case we we have for the annihilation channel 


|.M| 2 = 16 gN ( u +v' 2 ) 


{hpi? 


+ 


(hp2) 2 


m 2 x (kih) 


[ {Ml ~ t ) 2 (Mf - uf (M? - t)(M? - u) 


+ 2 w 


(2 


(fc IP]) 2 + ' rn l m X (pifcl) 2 + (pi fc 2) 2 - (^lh)(piP 2 ) + PPq rn \ 


(m 4 - ty 

(fclP2) 2 + 

(Mj - u y . 

+ 4 w , (u t + v ! )m q mx 


(Mi - t)(M | - u) 


+ 


+ 2 u'v' 
(fciPi) 


2 ™ 2 m 2 m 2 v rnl(p 1 p 2 ) 


rn q rn X , 


+ 


(Ml - ty (Ml - uy (Ml - t)(M} - tt) 
(fcip 2 ) (hpi) + (kiP 2 ) 


+ 


L( Ml-ty (Ml-uy 2(Ml-i)(Ml-u)\ 


+ 


(Zh - zhl 

cos 4 d w {m? z — s ) 2 


\{<?L + 4 ) (( fc iPi) 2 + (kiP2) 2 



c L c R m q ((P 1 P 2 ) - 2 m 2 x ) 


m 2 x (kik 2 )) 
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13 %u) 

2 cos 2 8 w {m 2 z — 5) 


(hpi) 2 + (hpiY 


(*?-«) (M\-u) 

+ ((Mf-T) + (ji^o) [(“'« - 

+ - t ,, ci)( 2 mj - (P1P2))] 

+ 2 m q mxw\cR — cl) 


(2v'cr — 2u cl) 


( 3 ) 


2(fcip 2 ) - (A 1 P 1 ) 2(fcjpi) - (kip 2 ) 


(M 2 d - u) 


(Ml - t) 


where mx is the neutralino mass, p\ and P2 are the incoming x four-momenta, k\ 
and ^2 are the outgoing fermion momenta, and (P1P2) denotes a four-momentum 
dot product. The symbols s, t, and u are the usual Mandelstam invariants, 
u f — a 2 + 6 2 , v f = a 2 + c 2 , and it/ = ea(b — c). 

The matrix element for the elastic scattering process XQ - ► XQ with momenta 
labeled p\ + k\ — > P2 + respectively can be found from eq. (3) by crossing 
&2 — * — &i, — ► ^2, and p2 — * — P2, so t — > u, u — > s, and 5 — > t. Likewise, 

the matrix element for the production process e + e” — ► XX (momenta labeled 
Pi + P2 — ♦ fci + £2) is found by <-+ p 2 and ^2 <-► pi. With these substitutions 
the labels m q and mx remain unchanged. 

To find the annihilation cross section we must integrate eq. (3) over the 
center-of-mass scattering angle. This is tedious due to the angular dependence 
in the propagator momenta, but we do not actually need this result because 
we are only interested in using the annihilation cross section to calculate relic 
abundances. For this purpose one needs a thermally averaged cross section, which 
is most easily obtained by expanding the relativistic cross section in powers of 
the relative velocity v . We do this before integrating and keep terms up to order 
v 2 . This method is also satisfactory for calculating neutralino annihilation in the 
Sun or galactic halo since again the collisions are non-relativistic. We find the 
annihilation cross section to be 
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tf'ann'y — ^ “ G^Cgmj^/3' j J/^ + V* 2 )(z 2 + (&1 + rj)i; 2 ) + 4tt/ 2 (l + (a2 + ^2)^ 2 ) 


where 


+ 4u/(u' + r')z(l + (ia: 2 + r 3 )?; 2 ) + 2uVz 2 (l + (a 4 + r 4 )u 2 ) 


+ (Z 2 3 - Z 2 4 ) 2 z 

+ (2, 2 3 - Z\, t )x'Y 


|( c i + c tf)(* 2 + “l” 2 ) + ^C£««2 2 (-l - atv 2 ) 

( v'cr - u'cl)(z 2 + (ai + z 2 rs)v 2 ) 


+ ( U ' C R ~ v’ci)z 2 (l + (a 4 + rs)t; 2 ) — 2 w'(ci — cr)z( 1 + (- x 2 + r 5 )v 2 ) 


v 2 0 l2 r ( -~{uc L + v'cr) + w'(c L - c R )z 



( 4 ) 


a 2 = -(2 — z 2 + x 2 ) 

1, ,, 

a 4 — — ( — 3 + x ) 

4 

V 

T\ = -(—4 + 2 2 + 4r — 3r^ 2 — rz 4 ) 

o 

r 2 = £(-5 + 2z 2 + 3r/3' 2 + 2 rz 2 /?' 2 ) 
r 3 = ^(-3 + 5r(3' 2 - 2/3' 2 ) 
r 4 = ^(-3 + 5 r/3' 2 ) 

rs = r -(- 3 - + rn 

where Gp is the Fermi constant, c q is a color factor, 3 for quarks, 1 for leptons, 
x> = Tn z/(( Tn 7 ~ -s ) 2 + ^m ^) 1 / 2 * s the P°i e factor, r z is the Z° width, 
z = Tn q /mx, 0' = (1 — z 2 ) 1 / 2 , and y' 2 = m^r/(M| + (3' 2 m 2 x ) is the squark mass 
suppression including the propagator momentum. The propagator momenta fac- 
tor r = m 2 x /(M? +P' 2 ) is usually small, as is a; 2 = \z 2 /(l- z 2 ). When m x — ♦ m, 
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however, z — ► 1, x — *■ oo, and the expansion breaks down. This occurs only very 
near mass thresholds, so we avoid evaluating eq. (4) right at threshold. The effect 
of ignoring propagator momenta can be found by taking the limits y' — + m\y jM q 
and r — » 0 in eq. (4). 

As a check, we can consider the pure photino (or pure higgsino) limit of 
eq. (4). A pure photino has Z \ 3 = Z\$ — a — 0, Z\\ — cos 0 W , Z 12 — sin^u, and 
b = c = e q sin 6 W . The annihilation cross section becomes 


0‘ann'y (photino) = )) 


8-k a 2 c q e*m^- (3' 

(Jlf? + m*./?' 2 )* 


z 2 + 


1 2 [-1 ^ 

-v 1 ; 

3 L 4 


, ^ 22 
+ -z x 


2 r + 2 r 2 - rz 2 + r 2 z 2 


3r 2 2 4 


( 5 ) 


5 3 2 19 

In the no propagator momenta limit we recover the well-known ’ 1 formula 

for the photino. We can also recover a higgsino by taking b = c = Z\\ — Z32 = 0, 
and Z\i — sin/3, Z\\ = cos/3. We perform the thermal average in the standard 
manner by replacing t; 2 with 6T jmx where T is the temperature, although this 
is not strictly correct to this order in x; 2 . 20 

9 

The above neutralino cross section, eq. (4), differs from previous work in 
several ways. First it includes the effect of propagator momenta which can be 
important if the neutralino is near in mass to the squark or slepton. Even for the 
pure photino, this effect can be important as is discussed in Sec. V and shown 
in Fig. 10. Second, several new terms appear which can be important or not 
depending on parameters. Third, we find opposite signs for the b 2 and c 2 terms, 
an important effect when Z° - squark interference occurs. 

To find the elastic scattering cross section we cross the matrix element as 
described above and integrate over the center-of-mass scattering angle. Since 
we are interested in the scattering of neutralinos from the halo of our galaxy, 
and the velocity of halo particles is known to be Vh a i Q ~ 10 -3 , we can take an 
extreme non-relativistic limit. This greatly simplifies the formulas. In particular, 
3 (m x + rn q ) 2 , t -► 2|pi| 2 (l - cos0*), u (m x - rn q ) 2 , and so if we can 
ignore the small mass of the quark, we have s % u and t ^ 0. Taking the limit 
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i 


O’el = 


8 G^m^rrij 2 
7t(t 71X + ”1g) 2 


(2a 2 + 6 2 + c 2 )a: 2 + -(c* - c £ )(Z 2 3 - Z 2 4 ) 


+ 4x 


ja 2 ( 6 - c) 2 |, 


( 6 ) 


where x 2 — — m 2 x ). We are actually more interested in scattering off 

nuclei than off quarks, and this will be discussed in the next section. Taking the 
pure photino limit we find 


<r c j(photino) = 


487ra 2 e^m^m g 2 
(m q + mx) 2 (A/| - m 2 x ) 2 ’ 


(7) 


which, apart from the newly included propagator correction to the scalar fermion 
mass, reproduces the well known formula. ’ A pure higgsino result can also be 
easily found. 

Finally, we wish to be able to include constraints on the parameter space 
from the ASP experiment. This experiment measured the cross section for e + e“ 
going to photon plus missing energy at the PEP storage ring at SLAC (yT = 29 
GeV). With their luminosity, the Standard Model background from e+e~ — > vu 7 
is predicted to be .03 pb, and from their limit of <r(e + e - — ♦ 7 + missing) < .06 pb, 
we get a limit on the production of neutralinos of a(e+e~ — ► XXl) < -03 pb . 24 In 
the soft photon limit valid here, the desired cross section can be calculated from 
cr(e + e“ — > XX) using the formula 26 


dcr(e + e 


lXX\a) 2a [(1 - |a;) 2 + \x 2 y 2 ] 


dxdy 


7 t x(l — y 2 ) 


cr(e 


+ e” 


XX] s) ( 8 ) 


where s is the Mandelstam variable, x = 2 E^jy/s is the dimensionless photon 
energy, y = cos 6 is the angle between the beam and the photon, and s = s(l — x). 
Performing the crossing described above, setting the electron mass to zero and 
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taking the limit Mgi = Mqr = Mg we find 


a ( 




e 




xj(6 4 + c 4 ) 


+ *5 X ' 2 ( Z 1 2 3 - ^14 )( C -R c ^ - Ci^ 2 ) + \x'\ Z h - ^l 2 4 ) 2 (ci + 4) , 

( 9 ) 

where x q = m\y/Mq and we have not included the propagator momenta here. 
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The pure photino limit of eq. (9) agrees with previous work. (See also Ref. 22). 


Ill, Effective Lagrangian 
and elastic scattering off nuclei 

The cross sections presented in Sec. II are for neutralino interactions with 
quarks or leptons, while for direct detection, or capture into the body of the 
Earth or Sun, the interaction between neutralinos and nuclei is needed. For this 
purpose, it is useful to re-derive some of the above using an effective Lagrangian 
approach. Since previous calculations have used this approach it is also useful 
for comparison purposes. 

We start from the supersymmetric Lagrangian interaction terms given in 
Refs. 7 and 8, and consider only the limit of heavy scalar fermions. We have 

— 4(7 2 ^ _ 

L eff = ~jj2~^ aPR + hP L)qq{ a PL + bP R )x 

M qL 

4 a 2 - 

- jj2~x{ cP R ~ aP L)qq{cPL - o-Pr)x ( 10 ) 

J qR 

2 

- ^2~( Z 13 “ Z u)^{ c lPl + C*Pji)gX7/i75X, 

where q is the quark field, and the other symbols were defined in Sec. II.. To get 
eq. (10) in a more useful form we perform Fierz transformations on the first two 
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terms. Using the fact that for Majorana fermions = 0 we find 


-£«// = ^2“|x7/i75X97 M (^' + 


^75)9 + 2ae(6sJ L - ci^)[xX99 + X75X?75?]|, 

( 11 ) 


where 


V = -;(<* + «i)(Z?3 - 2?4) + ^ t (i 2 - « 2 ) + * 2 S (« 2 - c 2 ) 

1 (12) 
= 2^ Cl ~ C *X Z 13 ~ Z h) - x \l(°? + i 2 ) - *J*(a 2 + c 2 ), 

and where = mw/M^i , etc.. Later we will include the effect of propagator 

momenta by substituting for etc.. For the limiting case of pure photino 
or pure higgsino eq. (11) agrees with the corresponding limiting cases of Ellis 
et al . For the general neutralino, however, there are several differences. First, 
there is a new term which is not of the form of an axial vector coupling. In the 
limit of equal left and right chiral squark masses, it is proportional to (b — c), 
the zino coupling, and to a, the higgsino coupling. This term may have been 
ignored previously because higgsinos, like Higgs bosons, couple proportional to 
mass, and for quarks or leptons this is usually quite small. However, for low 
energy elastic scattering off nuclei, Higgs bosons, 11 as well as higgsinos couple 
proportional to the nuclei mass, not the quark mass, and so this new scalar term 
can be important. Second, the sign of the b 2 and c 2 terms (or equivalently the 
(^13 ~ ^ 14 ) an d terms, since the overall sign is arbitrary) differ from Ellis et 
al.. Since in annihilation, for example, the a 2 term is usually small, the primary 
effect of this sign difference is negative Z°-squark interference, where before there 
was positive interference (and vice-versa). We shall show that this interference 
results in low direct detection rates for parts of parameter space. 

The effective Lagrangian, eq. (11) can be used to calculate the annihilation, 
elastic, and production cross sections in the standard manner. For the annihi- 
lation cross section there are several terms which appear in the complete cross 
section, eq.(4), which do not appear in the cross section calculated from the effec- 
tive Lagrangian. For the elastic cross section in the extreme non-relativistic limit 
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we find that the pseudo-scalar term does not contribute, and from the complete 
calculation, eq. (6), we also see that there is no interference between the axial 
vector and scalar terms. In fact, eq. (11) reproduces eq.(6) except for propagator 
momenta. 

We now use eq. (11) to find scattering off of nuclei. This is done by first finding 
^2 \(N\L e ff\N)\ 2 , where \N) is the nucleus wave function, which is assumed to 
be a sum of nucleon wave functions |n), which are in turn assumed to be sums of 
quark wave functions. The axial vector piece of the elastic cross section can be 
evaluated as in Goodman and Witten 1 and we find the elastic scattering cross 
section off a nucleus of mass mj\r to be 


2Am 2 x m 2 j S jG 2 p 

■K(mx + mx ) 2 



(13) 


where J is the total spin of the nucleus and the sums axe over the indicated 
quarks. The first term agrees with Ref. 1 in the photino limit (see also Ref. 19) 
and the second term is new and requires some explanation. In the above, we 
followed Goodman and Witten, 1 and Refs. 23 and 18 in defining 


A = 1 {1 + + 1 )-/(/ + 1 )V\J{3 + 1)]} (14) 

from the one particle nuclear shell model and the Lande formula, where / is the 
shell model angular momentum and s p is the proton (or neutron) spin. We also 
follow Refs. 16, 23 and 18 in defining (p\q^sq\p) = 2A qs qy where s q is the spin 
of quark q and A q measures the fraction of the proton spin carried by quark 5. 
Under some assumptions, the EMC group 1 gives Au = .746, Ad = —.508 and 
As = -.226, while the flavor SU(3) quark model predicts An = .97, Ad = —.28 
and As = 0. The EMC results are still controversial, so for most of our work 
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we use the flavor SU(3) values. In Sec. VII we show the effect of using the EMC 
values. Note that for simplicity we left out vector pieces in Eq. (11) which can be 
important if there is significant left and right chiral squark mixing. These terms 
have been discussed in Refs. 1 and 18 and are expected to be small. 

In deriving the second term of Eq. (13) we modified slightly a technique 

11 28 
described in Shifman et al. and used recently by Raby and West. For coherent 

scattering of a neutralino off a nucleus we need to find 


(N | 2 a(b - c)x 2 q qq\N) oc (N \ ^ Tzix 2 q d q m q qq\N ) , 

9 9 


(15) 


where (N\ is the nucleus state and the sum is over all the quarks, both valence 
and sea. Using the “heavy quark expansion” for the charm, bottom and top 
quarks Shifman et al. write Tn q qq ~ + O (^|) and by including the 

anomaly in the trace of the quark energy-momentum tensor 6^ they find 

Q s-y 

m N ^> N ^ N = {N\9^\N) ~ (N\G^G^\N) . (16) 


Physically, this last equation says that the mass of the nucleon (and therefore the 
nucleus) comes from the light quark anomaly. Since the light quarks in eq. (15) 
are very light, we can follow Shifman et al, in ignoring them and find 


(N\Y,^T 3L x 2 q qq\N) 

9 


4m jy \ v~> _ 2 , _ 2 / 2 Zh Z \ 3 \ 

27mw ) 3L q 1 27 mw \ s ^ n /^ cos /3 J ’ 


(17) 

where in the last step we made the simplifying assumption that all squarks have 
the same mass. Using Eq. (17), one finds Eq. (13) in a straightforward manner. 


The essence of the above derivation is that neutralinos can couple to the 
gluons in the nucleon via a loop involving heavy quarks and a squark. In the 
limit of very heavy squarks this becomes the same loop as Shifman et al. discuss 
for the Higgs exchange case. Since the higgsino, like the Higgs boson, couples 
proportional to the mass of the quark in the loop, the quark mass cancels out 
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and one is left with a coupling to the gluons which does not involve (to lowest 
order) the heavy quark mass. As mentioned above, the light quark glue gives 
rise to the mass of the nucleon, and so the final coupling is proportional not 
to the quark mass, but to the nucleon mass. For dark matter neutralinos from 
the galactic halo we can sum coherently over all the nucleons in the nuclei and 
arrive at a coupling proportional to the mass of the nucleus. Actually, when a 
heavy neutralino scatters off of a heavy nucleus some loss of coherence is possible 
and we take this effect into account in Sec. V. We do not claim that the above 
cross section, eq. (13) is exact, but it shows that “spin independent” cross sections 
exist for Majorana particles. Uncertainties include the extent to which the charm, 
bottom and top quarks contribute equally, the extent to which the strange quark 
contributes, the possibility of additional generations of quarks, and higher order 
contributions, both in the heavy quark and the heavy squark expansions. 

Finally note that we have included the effect of propagator momenta for the 
elastic cross section, eq. (13), by using x q , instead of x q in A! as derived in the 
complete elastic cross section. This can have a substantial effect if mx is near 
M q , as illustrated in Sec.V. 


IV. Relic Abundance and LSP constraints 

Using eq. (4) the present day mass density of neutralinos can be calculated. 
This “Lee- Weinberg” calculation has been the subject of many papers and we 
will use the method described in Refs. 30 and 31, which is an extension of the 
methods described in Ref. 32. This is an approximate analytic solution to the 
Boltzman equation n = —3 Hn — {<rv) ann (n 2 — n 2 ), which governs the number 
density n of particles as the universe cools. Here H is the Hubble parameter, 
no is the equilibrium number density of neutralinos and (&v) ann is the thermally 
averaged annihilation cross section derived in Sec. II. The basic idea is that 
when the temperature T of the radiation in the universe falls below the mass 
of the neutralino, the neutralino number density is suppressed by a Boltzman 
factor e~ E t T and falls rapidly. But since the interaction of neutralinos with 
ordinary matter is weak, there comes a time (denoted as freeze-out) after which 


18 


0 ‘"': •••••:. 

POOR 


L .c 


So 


Quality 


annihilation of neutralinos becomes improbable and the neutralinos present then 
are, for the most part, still present today. Our approximation uses an accurate 
method of finding the freeze out temperature, includes properly the effect of 
changing degrees of freedom and is accurate to better than 5% (as compared to a 
numerical integration). The reader is referred to Refs. 30 and 31 for more details. 

Figs. 2a,b,c show “scatter plots” of ft*h 2 versus the mass, mx, of the neu- 
tralino, where h = H/( 100 km sec -1 Mpc -1 ). These figures were made by 
selecting a value of Mq and tan/3 and considering a grid of points in the M, p. 
plane. For each set of the four parameters M, fx, tan/3, and we find mx and 
((Tv) , and then the relic abundance f Uh 2 . Each “x” in the figure represents 

a different supersymmetric model, and the collection of x’s indicate the range of 
relic abundances possible. There are several interesting features in these figures. 

First, notice that over a very wide range of supersymmetric parameters al- 
most no models give ft*h 2 < .01., (The lines across the figures indicate ft* = 1 for 
our preferred value h = |) the primary exception being at mx ~ mz/ 2, where 
annihilation is greatly enhanced due to the Z° pole. The values of M § and tan/3 
in Figs. 2a, 2c were selected because they give low values for ft*h 2 ; values of Mq 
of less than 50 GeV being inconsistent with results from the UA1 experiment at 
CERN. 13 The “wedge”’s on the left of the figures consist of models where the 
neutralino contains substantial higgsino content and move upward in ft*h 2 as 
tan/3 moves away from 1. The “lines” ’s through the Z° pole consist of models 
with a photino-like neutralino and also move upward in ft*h 2 , this time as the 
(assumed degenerate) scalar fermion masses increase. Recalling that the total 
luminous matter in the universe contributes ft ~ .01, the observed dynamical 
(dark) mass contributes ft ~ .1 — .3, and nucleosynthesis limits baryonic matter 
to .015 < fth 2 < .16, we see that if low energy supersymmetry exists and the neu- 
tralino is the LSP, it almost certainly is a significant component of the universe. 
This is true whether or not it is the main component of dark matter in galactic 
halos. Since mx is in the GeV range* neutralinos constitute cold dark matter 
and should cluster in galactic halos, so we consider it likely that galaxies contain 
as much or more mass in the form of neutralinos as in the form of baryons. This 
likelihood alone is enough to justify experimental attempts to detect neutralino 
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DM. It is also seen from Fig. 2 that a wide range of parameters predict a critical 
density of neutralinos (Q^h 2 = 1/4) and that much of parameter space is ruled 
out by requiring D* < 1. Since ft* = 1 arises naturally and since fi = 1 is 
desirable for several theoretical reasons, in this Report we will use the constraint 
x h 2 = 1/4 to reduce the parameter space which must be explored. As men- 
tioned in the introduction, galactic neutralinos may be detectable and many of 
the items discussed in this Report may be relevant even if < 1, especially 
since Fig. 2 suggests that Q* is likely to be greater than 4%. In fact, lowering 
requires an increase in coupling strengths causing neutralinos to interact more 
strongly in a detector. 

To begin to develop intuition about the neutralino parameter space we show 
in Figs. 3a,b,c,d contours of D-j/i 2 = 1/4 in the M,/x plane for several values of 
tan (3 and Mg. For typical values of tan (3 and Mg there is a “closed curve” in 
the middle of the figure and a “hyperbola” in the upper right hand corner. The 
contours are sometimes broken; small gaps resulting from the finite grid size in 
our contour finding program, and large gaps showing areas where no solution 
exists. No solution may exist when a neutralino mass eigenvalue crossing occurs 
and the resulting discontinuity in neutralino composition causes a discontinuity in 
relic abundance. For brevity, we will call a piece of broken contour a “strand” (as 
in spaghetti). These strands will be the building blocks of most of the figures in 
the remainder of this Report. We map the strands onto various other parameter 
spaces and use consistency and accelerator experiment results, which constrain 
the parameter space, to chop pieces off the strands. This is not, perhaps, the 
ideal way to present our results, but four (at least) dimensional parameter space 
is difficult to display. Note that Figs. 3a, b show the /x > 0 case while Figs. 3c, d 
show the /x < 0 case. 

As an important first example, consider the consistency constraint that the 
neutralino be the LSP , In a given model Mg is specified and so any values of M , 
fi, and tan/? which result in mx > M§ are ruled out. For example, all the models 
in Fig. 2a (with Mg = 50 GeV) to the right of mx = 50 GeV are inconsistent 
with the neutralino being the LSP and should not be considered. Throughout 
we discard regions of parameter space which do not satisfy mx < Mg. 
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Another important consistency constraint comes from the charginos, the su- 
persymmetric partners of the charged Higgs bosons. The neutralino parame- 
ters AT, and tan/3 determine uniquely the masses and couplings of the 
two charginos. Accelerator experiments constrain these masses and this will 
be discussed in Sec. VI, but here we only consider the consistency constraint 
rn X < rn x±' Figs. 4a, b show contours of mj — m^ ± for several values of tan/3 
and fi > 0. The areas between the two contours correspond to m^ ± < rax and 
are therefore inconsistent. Fig. 4a shows the mx = m x± contours on top of the 
= 1/4 contour for one value of Mg and tan 0. Note that much of one stand 
lies between the m x = r ^-x ± contours. This is a generic feature when /j, > 0 and 
cuts out considerable parameter space. (For /x < 0, we find m x < m i± almost 
everywhere and this constraint has little effect.) This loss of parameter space 
will show up later as breaks or gaps in the projected strands. Fig. 4b shows the 
m x — ± contours for several values of tan / 3 . 


V. Direct detection 


In this section we apply the elastic scattering cross section, eq. (13), to scat- 
tering from various elements and give the rates for neutralino interaction in a 
cryogenic detector. This new class of detectors plans to measure the small (order 
keV) energy deposited when a particle from the galactic halo hits a nucleus in 
the detector. Ionization detectors now operating with energy thresholds of order 
5 keV have already ruled out Dirac neutrinos with masses greater than 20 GeV 
as the major component of the galactic halo. 33 Many groups are developing new 
non-ionization cryogenic detectors which will operate at lower temperatures to 
reduce the background and lower the energy threshold. 

The rate of detection in a cryogenic detector is given in Ref. 3 (and verified 
by us) as 


R = 



( V )halo Phalo^el 
rn N m X 


(18) 


where cr el is the elastic cross section, eq. (13), (v) ~ 270 km/sec is the average 
dispersion velocity in the halo, r] v is a correction due to the motion of the Sun and 
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the Earth and rj t corrects for the energy threshold. Including only the motion 
of the Sun for simplicity, we find ~ 1.3. The factor r) c estimates the loss of 
coherence which occurs at high momentum transfer wh^mneutralinos start to 
interact with individual nucleons rather than the nucleus as a whole. Using the 
form factor from Ref. 3 


Vc-ot 


-2 


1 


. <i \ ~ 1 

, |2 , Hi. 

2 ^ y/Trexip' J 


-Y 


W (t^~) erf (th^) 

y/l + a 2 erfp' 


, (19) 


where 


a 2 ~ -83(^740 GeVg [(m^/200 GeV ) 1/3 + .06] 2 {{v) /270 km sec” 3 ), ( 20 ) 

{mx/m N + l) 2 i J 


and p f = ^\v 8 un! { v ) — 1 * For convenience we give a plot of r\ c as a function 
of ttijv and tti x in Fig. 5 . It is a substantial correction for heavy DM particles 
[rn x > 50 GeV) and for heavy nuclei (m N > 100 GeV). Note that for simplicity, 
we have taken a zero energy threshold throughout ( 77 $ = 1 ). The actual value of 
T] t depends upon detector design and is most important for light nuclei and light 
DM particles. 

In Fig. 6 we show the event rates in a mercury detector (natural abundance) 
plotted against the neutralino mass for several supersymmetric models. Fig. 6 a 
is the total rate, while Fig. 6 b shows the axial vector rate (without the new scalar 
term) and Fig. 6 c shows the scalar interaction rate alone. Figs. 6 a, b, c show the 
H> 0 case, while Figs. 6 d, e, f show the p < 0 case. We chose mercury not because 
it is an especially promising element, but because it emphasizes the possible 
importance of the new scalar term. It is heavy, and since A 2 J(J + 1 ) = 1/12 for 
the 17% of Hg which has a spin, it has very low rates for the pure axial vector 
coupling. Figure 7 is the same as Fig. 6 but for fluorine (7a, b), a light element 
which is 100 % a spin \ isotope with favorable shell parameters (A 2 J(J + 1) = 3/4) 
and for thallium ( 7 c), a heavy element with 100 % A 2 J(J + 1 ) = 3/4. For fluorine 
the axial vector term dominates almost everywhere while for thallium the scalar 
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and axial vector contributions are comparable. The lines in the figure correspond 
to “strands” in the M,/x plane for which 0*h 2 = 1/4 as discussed in Sec. IV. 
Various values of Mg and tan/3 are shown. All squark and slepton masses have 
been taken degenerate; the non-degenerate case will be discussed in Sec. VII. 

The rates shown here are easily seeded to other elements with the same one- 
particle-shell-model structure (Hg is a shell model “neutron”, while F and T1 are 
“proton” ’s). The axial vector rate is scaled by 

[AJ(J + l)]jy tun (ru Hg + rajr) 2 [n ( 21 ) 

[A J(J + 1)]h 6 ™Hg {™N + m x ) 2 /h 6 ’ 


where fx is the relative abundance of the relevant isotope, while the scalar term 
is scaled by 

m 3 N (m Hg m X ) (^c ) N (22) 

m^ g ( + m x) 2 ( 7 7c)Hg 

For germanium, which is modeled as a “neutron” for example, these factors are 
roughly 4 and .4 respectively for light neutralinos. For elements which are not 
describable by single-particle protons or neutrons the sum over A' Aq must be 
redone to get the axial vector rate. The scalar rate is still given by the formula 
above. 

Figures 6 and 7 show several important features. First note the wide range in 
rates possible at a given tyix which comes from considering the general neutralino 
rather than just the pure photino. (See Ref. 18 for another recent discussion of 
this point.) A large percentage of the bottoms of the strands in the loops in 
Fig. 3 represent almost pure photino states and these strands cluster in the large 
“pure photino” blobs seen at around 3 and 30 GeV. The rate for a pure photino 
depends only on Mq so there is one blob for Mq = 50 GeV and one for Mq = 125 
GeV. The variation in rate is perhaps over-emphasized with Hg where, as a shell 
model neutron, the axial vector rate is extremely small, due in part, to the quark 
model flavor SU(3) spin structure functions. (See Ref. 36 for the rates for Hg with 
EMC structure functions.) Thallium is also heavy, but is a shell model proton 
and has favorable shell model parameters and shows much less spread in the total 
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rate. However, the axial vector rate alone drops below 10~ 4 events/kg/day i n 
several places due to the cancellations described below. 

Next note the large dip in rate at mjf ~ m^/2 45 GeV. The Z° pole 

makes annihilation of 45 GeV neutralinos very efficient, so small values of the 
couplings are needed for = 1. These small couplings remain in the elastic cross 
section, where the Z° pole is not present, resulting in low rates. A neutralino 
with m x ^ m z /2 would be difficult to detect directly, though its suppression 
here implies an enhancement in decay or accelerator production. In Figure 
6b (axial vector term) there are also large dips in event rate at m x — 6 GeV 
and m x ~ 20 GeV. These are another generic feature and come about from 
negative Z -squark interference. (There is, in addition, cancellation among the 
terms in the sum over A f Aq .) The values m x ~ 6 and m x ~ 20 are not special, 
and as and tan/? are varied these cancellations occur for all values of m x . 
These low rates are mitigated by the contribution from the scalar term as seen 
in Figs. 6a and 7a. Even for fluorine, where the scalar term is small elsewhere, it 
dominates here and gives a minimum value for the elastic cross section. Another 
thing to note is that the pure photino blobs become lines as Mq is varied, and 
that they do not move from Fig. 6a to Fig. 6b showing that the scalar term does 
not contribute to pure photino elastic scattering. 

In order to get a feel for the four dimensional parameter space we show in 
Fig. 8 the effect of varying M i with tan/? fixed, and in Fig. 9 the effect of varying 
tan/? with Mq fixed. Keep in mind that reduction of the parameter space will 
result in further chopping of the strands when we include constraints from 
accelerator experiments in the next section. 

Finally, we note that the inclusion of propagator momenta has had a sig- 
nificant effect when m x neared Mq . The largest effect was in the annihilation 
cross section which decreased, thereby requiring an increased coupling strength 
for a given The increased coupling strength increased the event rate, which 
was further increased by the inclusion of propagator momenta.in the elastic cross 
section. For the general neutralino, the curves in Figs. 6 and 7 differ when prop- 
agator momenta are ignored, but the trend is hard to see. For the pure photino, 
however, the requirement = 1/4 fixes the squark mass as a function of m x , 
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which in turn fixes the event rate. In Fig. 10 we show the event rate in a flu- 
orine cryogenic detector both with and without propagator momenta included. 
At mx — 10 GeV, the rate is 30% higher with propagator momenta, at mx — 30 
GeV it is double, and by 90 GeV it is more than five times higher. If no new 
channels or particle thresholds appear, we note that the rate for pure photino 
detection reaches a minimum at about mx = 130 GeV and then starts to rise 
again as mx approaches Mq. 

VI. Accelerator constraints 

We have so far used cosmological predilection (f = 1) and consistency re- 
quirements (Mq,m^ ± > mx) to reduce the rather large parameter space in which 
neutralino dark mattes lives. In our selection of values of /x, M y tan/?, and Mq 
we have also implicitly used theoretical prejudice (tan (3 not too different from 
1, soft supersymmetry breaking parameters less than 1 TeV) and accelerator 
results (Mg > 50 GeV). In this section we will further restrict the parameter 
space by including the result from PETRA 14 that m^ ± > 23 GeV and the re- 
sult from ASP that a(e^e~ — > xXl) — *03 pb. The first experiment restricts 
parameter space because the same three parameters: /x, ilf, and tan/3 which de- 
termine the neutralino mass and couplings also determine the chargino mass and 
couplings. (The charginos are the supersymmetric partners of the tw T o charged 
Higgs bosons.) The ASP experiment is relevant because the e + e“ — ► XX cross 
section is related to the XX e+e - cross section. In fact, in general, accelerator 
experiments are important for particle DM detection since the DM particles may 
first be found there, and because negative results bear on the feasibility of direct 
detection. Since limits from experiments will continue to improve it is important 
to know what parts of parameter space are affected. 

Figure 11 shows contours of m% ± = 23 GeV in the plane for several 

values of tan/3. Fig. 11a shows /x > 0 and Fig. 11b shows fi < 0. The areas 
between the contours in Fig. 11a have m^ ± < 23 GeV and are ruled out, while 
the areas outside the contours are ruled out in Fig. lib. The effect this has on 
our spaghetti rate plots (Figs. 6a, d) is shown in Fig. 13, where the ruled out 
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areas are marked with “x”’s. For p > 0 m^ny of the models with mj < 10 
GeV have been eliminated. Pure photinos, however, are relatively unaffected by 
this constraint (Note the “blob” at 3 GeV is not “x”’ed out.) since they occur 
mostly near the M = 0 axis, and the affected areas between the contours in 
Fig. 11 a contain mostly models where the neutralino has substantial higgsino 
and zino components. We note that if the limit on the chargino mass was pushed 
to m^ ± > 40 GeV, most of the parameter space with m x < 20 GeV (except for 
pure photinos) would be ruled out. Conversely, if a chargino was discovered with 
m£ ± < 40 GeV we would have a good idea of the neutralino mass and couplings. 
Figure 13a shows the fi > 0 case, while Fig. 13b shows the < 0 case. For p < 0 
the chargino constraint is not nearly as important because much less area is ruled 
out. 

The ASP experiment , 15 performed at the PEP storage ring at SLAC, and sim- 
26 

ilar experiments, in principle actually produce neutralinos. These are neutrino 
counting experiments whose primary aim is to limit the cross section cr(e + e~ — + 
7 z/z/), where the photon is detected and the presence of the neutrinos is inferred 
from their missing energy. For neutralinos with mass below the ASP threshold 
(~ y/s/2 ~ 14.5 GeV), the corresponding process e+e“ — ► XXI can take place 
and the ASP limit of <r(e + e~ — > 7 +missing) < .06 pb (90% confidence level) can 
be used to rule out those areas of neutralino parameter space which would ex- 
ceed this. Actually, the neutrino production process must exist as a background 
(cr(e+e“ — ► i/ 1 / 7 ) ^ *03 pb) and the rate for production of exotic particles is 
then cr exot j c < .03 pb. The ASP group has used used this limit to produce re- 
strictions on pure photino parameter space, that is, the selectron mass, and find 

25 15 

> 62 GeV for mx = 0 . However, a Bayesian analysis of the same data 
results in the weaker bound <r exo ti c < .046 pb, giving > 58 GeV for mx = 0 . 
Both statistical methods have problems, and since a recent report combines 
Bayesian analysis and the results from other experiments to arrive at a number 
near 68 GeV for the selectron mass, we will use the stronger (and simpler) limit 
<r exotic < *03 pb here. 

We have integrated eq. ( 8 ) over the ASP acceptance 20 < d < 160 degrees, 
Emax > 12 GeV, and pp > .8 GeV/ c, and show in Fig. 12 contours of cr(e+e“ — ► 
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XXl) = -03 pb for fi > 0. The areas below the long contours and within the 
wedges at the right are ruled out at the 90% confidence level. The contours for 
the n < 0 case are very similar, the only differences being that the long contours 
extend straight through to the right edge, and that there are no “wedges’”s. In 
contrast to the PETRA chargino constraint, we see that for fj. > 0 this experiment 
affects mostly photino neutralinos and has little effect on the rest of parameter 
space. (See also the recent report of of Tata et al.. 22 ) The effect of the ASP 
limits on the event rate vs. mj plots of Fig. 6 is shown in Fig. 13 where the 
ruled out areas are marked with big boxes. For n > 0 the pure photino blob and 
a few other points are affected, while for fi < 0, virtually every neutralino with 
m X < 5 GeV is ruled out. We see that the PETRA and ASP constraints work in 
a complimentary fashion to rule out most very light (mx < 5 GeV) neutralinos. 
As new experiments are performed we expect these constraints to strengthen 
and eventually rule out (or discover) the light neutralino. Unfortunately, we see 
that current experiments give little information on heavier neutralinos and this 
situation is unlikely to change in the near future. 


VII. Model dependence and other uncertainties 

In this section we consider the effect of relaxing a few of the the many sim- 
plifying assumptions we made throughout the bulk of this Report. Specifically, 
we consider the effect of non-degenerate scalar fermion masses and the effect of 
using the EMC rather than the flavor SU(3) proton spin structure functions. 

The assumption of equal mass squarks and sleptons played an important 
role in producing reasonable cryogenic detector event rates. There are six left 
chiral squarks, six right chiral squarks and an equal number of sleptons, all of 
which, in principle, can have different masses. In addition, there can be off- 
diagonal terms in the mass matrices which result in mixing. These effects have 
been discussed previously, 1 ’ 18 so to limit the possibilities, we will here consider 
only splitting between the squarks and sleptons, assuming degeneracy among 
the squarks and among the sleptons themselves. The off-diagonal terms are 
expected to be small, and the effect of varying squark masses, while possibly 
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important, is very analogous to squark-slepton splitting. To get an idea of the 

.35 

range of models possible, we surveyed approximately 20 supersymmetric models 
which used renormalization group techniques to predict the scalar fermion mass 
splittings. These included “superstring inspired”, “heavy top supergravity”, and 
“light top supergravity” models among others. The range of splitting predicted 
varied substantially from model to model; however, all the models we surveyed 
were contained within two extreme cases: Mq = Mq and Mg — 3 Mj. We used 
Mq = throughout this Report, so we show here the case Mq = 3 Mf to 
hopefully bracket the effect of scalar fermion splitting. 

In Fig. 14 we show the total event rates in a mercury detector for Mq — 3 Mj 
and the same value of tan/3 displayed in Fig. 8. The values of Mq in Fig. 8 
correspond roughly to values of M- { in Fig. 14, with Mq in Fig. 14 being three 
time larger. The curves have moved and overall there is a reduction in rate by 
one to two orders of magnitude. As explained in the introduction this is to be 
expected. We see that substantial scalar fermion splitting, if it exists, would 
make neutralino dark matter detection more difficult. 

Recent measurements 16 of the spin structure of the proton can be interpreted 
as giving the surprising result that very little of the spin of the proton is carried 
by the quarks. While the measurement and interpretation are still controversial, 
the results, if true, have important implications for neutralino detection. In par- 
ticular, the matrix element (p\Q~f^75<l\p) = 2A qs 9 is found to differ substantially 
from the flavor SU(3) quark model expectations. The effect of this difference is 
illustrated in Fig. 15 where we show the analogues of Figs. 6b and 7a using EMC 
rather than flavor SU(3) values for the Ag’s. The rates are substantially changed 
whenever the axial vector term dominates. The scalar term is unaffected. Since 
the scalar term dominates the rate for mercury, in Fig. 15a we show the axial 
vector rate only (see Fig. 6b), For fluorine the axial vector term dominates so 
here we show the total rate (see Fig. 7a). We see that the magnitude and sign 
of the change caused by using EMC depends sensitively on both the neutralino 
composition and the nuclear shell model parameterization. Photino-like states 
(especially the pure photino blobs) seem to be strongly affected, the rate be- 
ing either strongly enhanced for shell model neutrons like Hg, or suppressed for 
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shell model protons like F. Note that the positions of the Z -squark interference 
cancellations have moved and that for shell model neutrons there is a general low- 
ering of the rate. The effect of the EMC structure functions on direct detection 
of neutralinos has also been discussed recently in Ref. 18 . This uncertainty in 
spin structure functions will have to be resolved before accurate determinations 
of neutralino detection rates can be made. 

Finally we mention some other uncertainties and limitations bearing on the 
results of this Report. First, we note we have ignored the top quark and the 
Higgs bosons throughout. These particles almost certainly exist in some form, 
but their masses and properties are unknown. The cross sections and some of our 
conclusions would change for some values of the top mass. Ellis et aL in fact, 
included a 30 GeV top quark in their original calculations, and this explains some 
of our differences. The Higgs bosons (there are five in the minimal supersymmet- 
ric models) are even more problematic since both masses and mixing parameters 
are unknown and again, their inclusion could make important changes in our 
results. Next, we have throughout only considered neutralinos less massive than 
the Z° boson. Heavier neutralinos could exist, but new annihilation channels 
open up and the problem of ignoring the top quark and Higgs bosons is exac- 
erbated. Certainly, the trend is that detection becomes more difficult as the 
neutralino mass increases, but a very heavy LSP pushes the masses of the other 
supersymmetric particles even higher and makes a supersymmetric solution to 
the hierarchy problem more unlikely. Another limitation is that we considered 
only the minimal supersymmetric extension of the Standard Model. Although 
this is the most well studied class of models, many other models have been pro- 
posed, and the addition of new particles and couplings can change our predictions 
drastically. In fact, there may even be particle dark matter whose only interac- 
tion with ordinary matter is gravitational. Finally, and most importantly, we 
have throughout made the assumptions that low-energy supersymmetry exists, 
and that the neutralino is the LSP. 


VIII. Conclusions 

In this Report we covered a wide range of topics having to do with the de- 
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tection of neutralino dark matter. We started by calculating the neutralino an- 
nihilation, elastic scattering, and production cross sections and then used these 
to find the relic abundance of neutralinos and the rate of interaction of galactic 
neutralinos with a cryogenic detector. New features included additional terms in 
the cross sections, the inclusion of propagator momenta, and the sign of inter- 
ference terms. We emphasized that over most of the supersymmetric par am eter 
space a near critical density of neutralinos exists, and that therefore they are 
likely to constitute a significant fraction (ft* > .04) of the galactic halo. This is 
true whether or not they are, in fact, the dark matter. With regards to direct 
detection, a new term in the elastic cross section coming from heavy quark loops 
gives neutralinos a “spin independent” interaction which can be important for 
heavy materials, and reduces the Z°-squark negative interference in almost all 
cases. Pure photinos and pure higgsinos are not affected by the new scalar term. 

We then considered the complementary nature of direct detection DM searches 
and accelerator experiments, showing how results from the PETRA, ASP, and 
UA1 experiments rule out large areas of parameter space. These constraints in 
turn affect the detectability of neutralinos. Finally, we considered the effect of 
non-degenerate squarks and sleptons on our results, showing a reduction in direct 
detection rates for the most extreme cases of one or two orders of magnitude, 
and also the effect of the new EMC proton spin structure functions on the rates, 
showing again a substantial change. 

Overall, we found that the neutralino makes a superb particle dark matter 
candidate, and that the pure photino and pure higgsino special cases usually 
considered do not represent well the breadth of possibilities. These particles have 
a chance of being detectable, either directly or in accelerator experiments and we 
encourage our experimental colleagues to make strong efforts in this direction. 
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Appendix A 

In this appendix we list the complete matrix element, including propagator 
momenta and non-degenerate squark (or slepton) masses. We also show the 
initial matrix elements corresponding to the Feynman diagrams in Fig. 1. 

For the annihilation process (xx 9?) labeled by the four-momenta (pi + 
P2 + k 2 ) the five graphs of Fig. 1 have matrix elements 


M a 

M b 

M c 




- 0 
~2g 2 

~2g 2 

\M 2 l - u) 


u(k\){atjP L + bPR)u(pi)v(p 2 )(aeiP R + bPi)v(k 2 ) 
u{h)(-aP R + ejcP L )u(p 1 )v(p 2 )(-aP L + ticP R )v(k 2 ) 
v(p 2 )(aP L + ejbP R )v(k 1 )u(k 2 )(aP R + tjbP L )u{pi) (Al) 


Mj = 


-2g 7 


("?* - u) 


v(p 2 )(-ejaP R + cPi)v(ki)u(k 2 )(-eiaP L + cP^)u(pi) 


Mz = 2 ~ co S 2 el\m 2 z- !) V ^ P2 ^^ 5U ^ Pl W 1 * 1 h^( c L P L +c R P R )v(k 2 ), 


where e, is the sign of the neutralino mass eigenvalue and the other symbols are 
defined in Sec. II, 

The total matrix element is M — M z + M c + — M a — The total 

matrix element squared is the sum of the following terms. 


\M C 


16 g 4 


(a 1 + b 2 ) 2 (kipi)(k 2 p 2 ) + Aa 2 b 2 m 2 mx 


(Mj L ~ *) 2 

+ 2 m q mxab(a 2 + b 2 )e((kipi) + ( k 2 p 2 )) 

t y ( fl2 + c2 ) 2 (^lPl)(^2P2) + 4o 2 c 2 m, 2 my 5 
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\M C \ : 


— 2 m q mxac(a 2 + c 2 )e ((&ipi )(fc 2 p 2 )) 
16 5 4 


(A2) 


(a 2 + fc 2 ) 2 (fcip 2 )(fc 2 pi) + 4 a 2 b 2 m 2 x m q 


\M d \ 2 


(m? r - uy [ 

+ 2 m q mxab(a 2 + & 2 )e((&ip 2 ) + (fopi)) 

16 5 4 


(a 2 + c 2 ) 2 (fcip 2 )(fc 2 pi) + 4a 2 c 2 m q 2 m 2 x 


(MIr-u) 

— 2 m q mxac(a 2 + c 2 )e((fcip 2 ) + (fopi)) 


\M Z \ 2 = 


16g 4 (Z 2 3 — ^ 2 4 ) 2 
cos 4 6 w (m 2 z — s ) 2 [4 


t (4 + c p)(( fc iP2)(fc2Pi) + (*lPi)(fc2P2) - rn 2 x {kyk 2 )) 


+ ^c i c jR 77 i g 2 ((pip 2 ) - 2 m|) 


2 f 2 eA 4 a .M£ 


32 p 4 


a(c — 6)(fcipi) + m q m x e(bc — a 2 )] 


(M? l - t)(Mf s - () 

x [a(c — b)(k 2 p 2 ) + Tn q m x e(bc — a 2 )] 


2ReM a M* = 


16 g 4 


-«){ 


2 a 2 6 2 


(pifc 2 )(p2*;i) 


( M q L - i){M q L 

+ (fciPi)(^2P2) - {k\k 2 ){pip 2 ) + Tn q 2 m 2 x + m g 2 (pip 2 ) 
+ (a 4 6 4 )m^(A:iA: 2 ) 

+ m, q mxabe(a 2 + b 2 )[(p 2 &i) + (pifc 2 ) + (pi^i) + (p 2 fc 2 )] | 


2 i 2 eA 4 a A 4 ^ = 


16g 4 [ 


2 a 2 be 


(M qL - t)(M qR 

- {k\k 2 ){p\P7) + m q 2 m 2 x 


{P\k2){P2k\) + (fclPl)(fc 2 P 2 ) 


+ a 2 (b 2 + c 2 )m 2 x {k\k 2 ) + m g 2 (a 4 + fc 2 c 2 )(pip 2 ) 

+ m q mxae[b(a 2 + c 2 )((p 2 fci) + (pik 2 )) - c(a 2 + & 2 )((pi&i) + ( p 2 k 2 ))] j 
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2ReM b M* 


2ReM b M d 


2 ReM c M* d 


2ReM a M*. 


2ReM b M* 


16g 4 


—2 a 2 bc 


{p\k 2 ){p 2 ki) + {kipi)(k 2 p2) 


( MqR - t)(MqL - U ) 

- ( kik 2 )(piP 2 ) + m q 2 m x 
+ a 2 (b 2 + c 2 )rn 2 x {k\k 2 ) + m q 2 (a 4 + 6 2 c 2 )(pip 2 ) 

+ m q mxae[-c(a 2 -f b 2 ){(p 2 ki) + ( p\k 2 )) + 6(a 2 + c 2 )((pi 
16flr 4 



2aV 


(pik 2 )(p 2 ki) + (fcipi )(^2P2 ) 


{^qR ~ i){M q R ~ U ) 

- (kik 2 ){pip 2 ) + m q 2 m 2 x + m, 2 (pip 2 ) 

+ (a 4 + c A )m 2 x {kyk 2 ) 

— m ? mjace( a 2 + c 2 ) [(p 2 ) + (pi& 2 ) + (pi&i) + (P2&2)] | 


32p< 


a(c — &)(&ip 2 ) m q m x e(bc — a 2 )] 


( M 2 L -u)(M 2 R -u) 

x [a(c — b){k 2 p\) + m q m x e(bc — a 2 )] 
8g\Z^ - Z 2 14 ) 


2 (6 ex - a 1 c R )(kipi)(k 2 p 2 ) 


cos 2 ~ s ){M~i — t) 

+ 2 {a 2 ci - b 2 CR)m q 2 m 2 x - ( a 2 C£ - b 2 CR)m q 2 (pip 2 ) 
- ( b 2 c L - a 2 CR)m;r(fci& 2 ) 

+ abtm q m x {cR - c L )[(pik 2 ) + (p 2 fci) - 2(fcipi) - 
%9 A {Z 13 ~ %u) 



2ReM c M* z 


cos* «„(m| - , ){M > r > 2 <“ c <=*)(*>«)(*»«) 

+ 2(c 2 cx - a 2 CR)m q 2 m 2 x - ( c 2 C£ - a 2 c^)m 9 2 (pip 2 ) 

- (a 2 Ci - c 2 c^)m^(A;]fc 2 ) 

+ acem q m x {cR - c L )[{-pik 2 ) - (p 2 fci) + 2 (fc!pa) + 

W (Ik- 7)(Affc _ „) { 2 (“ 2c * - V'lVMirtbTi) 

+ 2 ( 6 2 cr - a 2 c L )m q 2 rn 2 x - [b 2 cr - a 2 CL)m q 2 {p\p 2 ) 
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2 ReM d M* z 


- ( a 2 c R - b 2 CL)m 2 x (kik 2 ) 

+ abem q mx(cR — ci)[2(pik 2 ) + 2(p 2 ki) — (fcipi) — (& 2 P 2 )] 


&9 4 {Z 13 ^14 ) 


2(c 2 cr — a 2 cx,)(fcip2)(fc2Pi) 


cos 2 6 w (m 2 z — s)(M? r — u) 

+ 2(o 2 cr - c 2 c L )m q 2 m 2 x - ( a 2 CR - c 2 c L )m q 2 (pip 2 ) 

- (c 2 cjj - a 2 Ci)m^-(A:ifc2) 

+ actm q m x (cR — cl) [-2(pifc2) — 2(p2&i) + (^lPi) + (&2P2)] |> 


where all symbols were defined in Sec. II. 
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Figure Captions 

1. Feynman diagrams contributing to neutralino (x) interaction with quarks 
or leptons ( q ). The first diagram shows Z° exchange while the last four 
show left and right chiral squark (or slepton) exchange. 

2. Relic abundances of neutralinos vs. the neutralino mass. Each “x” shows a 
different supersymmetric model. The horizontal line denotes critical density 
(ft = 1) for h = Supersymmetric models with parameters in the wide 
range 0 < M < 1 TeV and 0 < // < 1 TeV are shown for several values of 
squark mass M$ and tan /?. 

3. Contours of ft^h 2 = | in the M, n plane for various values of Mg and 
tan/?. Figs. 3a and 3b show the fi > 0 case while Figs. 3c and 3d show 
/i < 0. Solid lines indicate tan/? = dashed lines indicate tan/3 = 2 and 
dot-dashed lines indicate tan fi = j. 

4. Consistency constraint: contours of m x = m i± in the M, fi plane for 
various values of Mg and tan/?, where m i± is the mass of the lightest 
chargino. The areas between the contours are inconsistent with the neu- 
tralino as LSP. Fig. 4a shows the ft^h 2 = | contour (solid lines) being cut 
by the m x = contour (dashed lines) for Mg = 50 GeV and tan/3 = 
Fig. 4b shows contours for Mg = 100 GeV and several values of tan/?. See 
Fig. 3 caption for the code for lines. 

5. Coherence loss factor r\ c as a function of neutralino and nuclei mass. The 
various lines are labeled on the right by the nuclei mass in GeV. 

6. Event rates in a cryogenic detector made of mercury as a function of neu- 
tralino mass. All points correspond to ft^h 2 = 1 for tan/? = .25, .75, or 2 
and M$ = 50 or 125 GeV . Lines in Figs. 6a, d are labeled by tan fi. The solid 
lines indicate M$ = 50 GeV while the dashed lines indicate M$ = 125 GeV. 
Figs. 6a,b,c show the p. > 0 case while Figs. 6d,e,f show the p, < 0 case. 
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Figs. 6a, d show the total rate, Figs. 6b, e show the contribution from the 
axial vector term alone and Figs. 6c, f show the new scalar term contribution 
alone. 

7. Same as Fig. 6 for detectors made of fluorine and thallium. Only the fi > 0 
cases are shown. 

8. Effect of varying squark mass: Event rates in a thallium detector as a 
function of neutralino mass for tan/? = | and Mq = 50, 60, 75, 100, 125, 
150, and 200 GeV. The lines and “pure photino blobs” are labeled with the 
squark mass in GeV. For fi > 0 (Fig. 8a) no blobs were found for Mq ^ 125 
geV, and for /x < 0 (Fig. 8b) no blob was found for Mq = 200 GeV. 

9. Effect of varying tan/?: Event rates in a thallium detector as a function of 
neutralino mass for Mq = 50 and tan/? = .25, .33, .75, 1, 3. The lines are 
labeled by tan/?. 

10. Effect of propagator momenta for the pure photino. For each value of 

photino mass ttixi Mq I s f° un d so that ^ and the corresponding 

event rate in a fluorine cryogenic detector displayed. Curves show the effect 
of propagator momenta both in the annihilation cross section and in the 
elastic scattering cross section. Flavor SU(3) spin structure functions were 
used. 

11. PETRA constraint. Contours of m i± = 23 GeV in the M, /i plane are 
shown for several values of tan/?. Fig. 11a shows the fi > 0 case where the 
areas between the contours are ruled out due to Tn-£ ± < 23 GeV. Fig. lib 
shows the \jl < 0 case where the areas outside the contours are ruled out. 

12. ASP constraint. Contours of &ASP = -03 pb (see text) in the M , /i plane 
are shown for Mq — 50 and several values of tan/?. The code for the lines 
is the same as in Fig. 3. The areas below the long contours and inside the 
“wedges” at the far right are ruled out. Only the fj, > 0 case is shown; the 
fi < 0 case is very similar except the horizontal contours continue straight 
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across to /x = 1 TeV and there are no “wedges”. There are no solutions for 
Mq — 125 GeV and so no constraints exist for that case. 

13. Effect of PETRA and ASP constraints. Event rates in a mercury detector 
as a function of the neutralino mass (same as Figs. 6a, d) with areas ruled 
out by the PETRA and ASP constraints marked. Fig. 13a shows the fi > 0 
case while Fig. 13b shows the /x < 0 case. Areas ruled out by PETRA 
are marked with “x” ’s and areas ruled out by ASP axe marked with large 
boxes. 

14. Effect of non- degenerate squarks and sleptons. Event rates in a thallium 
detector as a function of the neutralino mass with — 3Hf^ instead of 
Mq = Afj- (see Fig. 8). Lines and blobs are labeled by the value of the 
squark mass in GeV. In comparing with Fig. 8 the relevant mass is actually 
the selectron mass; one third of the labeled mass. Fig. 14a shows the /x > 0 
case while Fig. 14b shows the /x < 0 case. 

15. Effect of EMC structure functions. Same as Fig. 6b (15a) and Fig. 7a (15b) 
with EMC rather than flavor SU(3) spin dependent structure functions. 
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